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Abstract. We study, at T = 0, the anomalies in the magnetization curve of 
the 5 = 1 two-leg ladder with frustrated interactions. We focus mainly on the 
existence of the M = M B /2 plateau, where M B is the saturation magnetization. 
We reports the results by degenerate perturbation theory and and the density 
matrix renormalization group, which lead to the consistent conclusion with each 
other. We also touch on the M = M B /A and M = (3/4)M s plateaux and cusps. 
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Anomalies in the magnetization process, such as plateau, cusp and jump, in low- 
dimensional magnets have been attracting increasing attention in these years. In this 
paper, we investigate the effect of the frustrated interactions on the plateaux and cusps 
in the magnetization curve of the 5 = 1 two- leg ladder. Our Hamiltonian, sketched 
in Figure 1, is described by 
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where Sjj is the S = 1 operator at the jth site of the Zth ladder (I = 1, 2), H denotes 
the magnetic field along the z direction, and all the couplings are supposed to be 
antiferromagnetic unless otherwise noticed. 
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Figure 1. Sketch of the model Hamiltonian iTi . 

We focus mainly on the effect of the frustrated interaction J 2 on the M — M s /2 
plateau, where M is the magnetization and M s is the saturation magnetization, 
although we touch on other topics. In §2, we investigate the simple two-leg ladder 
with no frustration. The effect of the frustrated interaction J 2 on the M = M s /2 
plateau is discussed in §3. The last section §4 is devoted to concluding remarks. 

2. M s /2 plateau of simple 5=1 two-leg ladder 

In case of the 5 = 1/2 simple ladder, the rung interaction is always relevant. In other 
words, infinitesimally small rung interaction brings about the spin gap at M — 0. The 
situation is quite different for the M = M s /2 plateau of the present 5 = 1 model. 
When J — * 00, the M — M s /2 plateau obviously exists, because the problem is 
reduced to the two-spin problem. On the other hand, there will be no M = M s /2 
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plateau in the Jo — > — oo limit, because the ladder is essentially the single chain of 
5 = 2 spins formed by the rung spin pair. Thus there exists the critical value of 
Jo/Ji. This quantum phase transition is thought to be of the Berezinskii-Kosterlitz- 
Thouless (BKT) type [1,2]. We note that this M s /2 plateau state is unique (not 
degenerate) from the necessary condition for the plateau [3] . It is interesting whether 
the critical point lies in the antiferromagnetic side (i. e., Jo > 0) or the ferromagnetic 
side (Jo < 0), when Ji > is fixed. Infinitesimally small Jo yields the M = M s /2 
plateau if Jq < 0, while it does not if Jq CI ^ > 0. 

To estimate the above-mentioned critical point analytically, we employ the 
degenerate perturbation theory (DPT) [4]. Hereafter we set Jo = 1 (energy unit) 
for convenience and due to the fact that the critical point lies in the antiferromagnetic 
side, as will be seen later. Let us begin with the strong rung coupling limit Ji <C Jo. 
In this limit, around M = M s /4, we can take only two states for the rung states, 
neglecting other 7 states: the lowest state with S'rimg*' 1 = and that with S'mng*'' = 1. 
We can express these states by the T z — 1/2 and T z = —1/2 states, respectively, 
of the psuedo spin T. The lowest order perturbation with respect to Ji leads to the 
effective Hamiltonian 

^eff - £ K t ;V/ f/ ;\; + T]T] +1 ) + J* S T*T* +1 - H cS T*} (2) 
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J3 = ^> = H eS = H-l-±. (3) 



The M = 0,M s /4 and M s /2 states of the original S system correspond to the 
M (T *> = -A/j T) ,0 and Af s (T) states, respectively, where M (T *> (M S (T) ) denotes the 
magnetization (saturation magnetization) of the T system. It is easy to obtain the 
field corresponding to = by considering the one-spin-down spectrum of 

Ti-eS as 

rr(!) _ i , 11J 1 (A\ 
H M s /2 - 1 + -3— (4 ) 

This , 2 gives the lower edge of the M = M s /2 plateau. 

Similar DPT can be developed around M — (3/4) M e , resulting in 

<; /2 =2-Ji, (5) 

where H y M ' , 2 gives the upper edge of the M — M s /2 plateau. Thus, the critical value 

Ji where the M = M s /2 plateau vanishes can be estimated from = H^^ 2 , 
resulting in 

4 CT) = A = 0.214. (6) 
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We remark that there is no M = M s /2 plateau for the so-called isotropic case J\ = 1. 

We have calculated the magnetization curves (Figure 2) by use of the density 
matrix renormalization group (DMRG) method. The DMRG result is consistent 

(cr) 

with the LS result Jj ' — 0.491, if we consider the pathological nature of the BKT 
transition. We can clearly see the M s /2 plateau J\ = 0.3, while we cannot for J x = 0.5. 
It is difficult to judge for the Ji = 0.4 case. This situation is qualitatively consistent 
with the DPT results, although the critical value is slightly larger than the DPT 
prediction ©. This is reasonable because the plateau region will narrowly extend like 




Figure 2. Magnetization curves of the simple ladder (J2 = J3 = 0) when 
Ji = 0.3, 0.4 and 0.5 for 40 spins obtained by the DMRG. The M s /2 plateau 
clearly exists for Ji = 0.3, while it does not for Ji = 0.5. It is difficult to judge 
for the Ji = 0.4 case. 



a beak of a bird on the J\ — H plane, as was seen in the diamond type spin chain 
case [5]. 

The level spectroscopy (LS) method [6-8] is very powerful in finding the quantum 
critical point of the BKT and Gaussian types. We have also performed the LS, the 
details of which will be published elsewhere. Our conclusion is 

j[ a) = 0.491. (7) 

which is consistent with the DMRG result. We have also performed the non-Abelian 
bosonization approach, finding a qualitatively consistent conclusion with those by the 
above methods. Details of the LS method and the non-Abelian bosonization approach 
will be published elsewhere. 



3. Effect of frustrated interactions 



Let us consider the effect of J 2 interactions on the M s /2 plateau problem. The DPT 
in the previous section can be easily extended to the case with J 2 . We obtain 

< ) / 2 = l + 1 ^-^ < /2 = 2-J 1+ 3J 2 , (8) 
yielding the critical line 

J2 = Ji ~ y A (9) 




Figure 3. Magnetization curve for 60 spins when J\ = 0.5 and J2 = 0.4. 

The magnetization curve for the J\ = 0.5, J2 = 0.4 case calculated by the DMRG 
method is shown in Figure There exists the M s /2 plateau which is not seen in the 
J% = 0.5, J2 = case (no frustration case), as shown in Figure[21 Furthermore we can 
see the M s /4 and (3/4)M s plateaux, which is attributed to the Neel ordering of the 
T system, which was first pointed out by our group [9]. The Neel ordering condition 
is known from 7i ff as J^g < J* s , from which we obtain 

■h > ^Ji = 0.684J 1; (10) 

for the M s /4 plateau. This condition is satisfied for the parameter set of Figure |21 
The Neel ordering condition for the (3/4)M s plateau is J2 > O.6J1, which explains 
the fact that the width of the M = (3/4)M s plateau is wider than that of the M s /A 
plateau. Both of the M s /A and (3/4)M s plateaux require the spontaneous breaking of 
the translation symmetry as is known from the necessary condition for the plateau [3] . 
We have also employed the LS for this case, finding consistent results with that of the 
DMRG. 



4. Concluding remarks 

Here we shortly touch on the effect of J3. In the framework of the DPT, the J3 
interaction brings about the frustrated next-nearest-neighbor interaction between Tj 
and Tj + 2, while the J2 interaction does not. Thus, for sufficiently large J3, the 
magnetization curve has cusps as found in Figure 5. These cusps are due to the 
mechanism proposed by Okunishi, Hieida and Akutsu [10]. The condition for the 
existence of the M s /4 plateau is J 3 > 0.31Ji from the DPT [9]. The width of the 
plateau may be too narrow to be observed clearly in Figure 0| 

We have discussed the plateaux and cusps in the magnetization curve of the 
5 = 1 frustrated two-leg ladder by use of the DPT as well as the DMRG. We have 
also used the non-Abelian bosonization approach and the LS method, although we 
did not enter into their details. The results obtained by these methods are consistent 
with each others. We remark that the M s /A plateau of the present mechanism is 
possibly related to that observed in organic 5=1 spin ladder 3,3',5,5'-tetrakis(./V- 
£eri-butylaminxyl)biphenyl (abbreviated as BIP-TENO) [11]. 




Figure 4. Magnetization curve for 40 spins when Ji = 0.3 J2 = and J3 = 0.2. 
We can see cusps near M = 0.15M S and M = 0.85M S . 
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